o 



X 



UMN-TH-2832/10 

FTPI-MINN-10/03 

January 2010 



The violent Universe: the Big Bang* 



K. A. Olive 

William I. Fine Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455, USA 

Abstract 

Four lectures on Big Bang cosmology, including microwave background radi- 
ation, Big Bang nucleosynthesis, dark matter, inflation, and baryogenesis. 

The Big Bang theory provides a detailed description of the history and evolution of the Universe. 
Direct experimental and observational evidence allows us to probe back to the first second after the 
initial state (bang) when the temperatures were of order 1 MeV and the light elements were created. Our 
understanding of the Standard Model of electroweak interactions allows us to push the description of the 
early universe back to about 10^^'^ s after the bang when we expect that the electroweak symmetry was 
j^. restored. Indeed, it is possible to discuss events back to the Planck time (10~^^ s after the bang) albeit in 

a very model dependent way. 

In these four lectures, I hope to give an overview of modern cosmology with an emphasis on 
^ I particle interactions in cosmology. After a description of the standard FLRW model (including the 

microwave background radiation) in Lecture 1, 1 will cover both inflation and baryogenesis in Lecture 2. 
' (— I ' , Lecture 3 will focus on Big Bang Nucleosynthesis (BBN) and Lecture 4 on dark matter. 

0-<\ 1 Lecture 1: Standard Cosmology 

^ . 1.1 The FLRW metric and its consequences 

The standard Big Bang model assumes homogeneity and isotropy. As a result, one can construct the 
►^ , space-time metric by embedding a maximally symmetric three dimensional space in a four dimensional 

\f^ \ space-time (see, e.g., Ref. HI). The most general form for a metric of this type is the Friedmann- 

'^ ■ Lemaitre-Robertson-Walker metric which in co-moving coordinates is given by 



in- ds^ = dr - R\t) 

o 
o 



(1 - kr^) 



(1) 



where R{t) is the cosmological scale factor and k the three-space curvature constant {k = 0, +1, —1 
for a spatially flat, closed or open universe), k and R are the only two quantities in the metric which 
distinguish it from flat Minkowski space. It is also common to assume the perfect fluid form for the 
energy-momentum tensor 

P_ ; T'^^ = vg^" + {V + p)u''u'', (2) 

where g^,y is the space-time metric described by Eq. ([Til, p is the isotropic pressure, p is the energy density 
and u^ = (1,0,0,0) is the velocity vector for the isotropic fluid. The (00) component of Einstein's 
equation 

Rfiu - -^g^iuR - ^9fMu = SttGnT^u (3) 



yields the Friedmann equation 




^,^,.| ^8^_|.^A_ ,,, 
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and the {ii) components give 



(5) 



3 3 

where A is the cosmological constant. In addition, from T^'^-i, = 0, we obtain 

p=-3H{p + p). (6) 

Note that of these last three equations, only two are actually independent. These equations form the basis 
of the standard Big Bang model. 

At early times (t < 10^ yrs) the Universe is thought to have been dominated by radiation so that 
the equation of state can be given by p = p/3. If we neglect the contributions to H from k and A (this is 
always a good approximation for small enough R) then we find that 

R{t) ~ t^/2 (-7) 

and p ~ R~^ so that t ~ {3/32ttGnp)^^'^- Similarly for a matter or dust dominated universe with p = 0, 

R{t) ~ t^/s (-g) 

and p ~ R^^. The Universe makes the transition between radiation and matter domination when prad = 
Pmatter or whcn T ~ fcw X 10^ K. In a vacuum or A dominated universe (that we are approaching 
today) 

i2(t)~e\/^*. (9) 

More general solutions for the behaviour of the scale factor can easily be found by defining a 

quantity Q: 

3k 
Q = J^- SttGnP. (10) 

If we further assume p = (7 — l)p (with 7 between 1 and 2), we have that p ~ R^^'^ and we can write 

1/2 



R 

R 



A-Q 



(11) 



implying that Q < A. For all choices of k, the function Q — )• —00 as i? — )• 0. When k = +1, it is 
easy to see that Q has a maximum (with Q > 0) and then tends to as i? — )• 00. When A; = or —1, 
Q monotonically increases towards as i? — ;■ 00. In Fig. \T\ the function Q is plotted qualitatively as a 
function of the scale factor. 

Let us first consider the more interesting case of A; = +1. For a cosmological constant A > Qmax, 
Q < Aior all R. In this case, two distinct solutions are possible: the universe may expand forever from 
a singularity at i? = to i? = 00, or by choosing the lower sign in Eq. (ITTI ). the universe may collapse 
from infinity to a singularity. 

There are also two possible solutions for < A < Qmax, the case depicted in Fig.[T] The universe 
may again start at a singularity at i? = and expand to the point where Q = A and then recoUapse. 
Alternatively, the Universe may start from infinity and collapse to the point where Q = A (at larger 
R), bounce back and expand to infinity. When A = Qmax there are a total of five solutions. At the 
value of R such that Q = Qmax, we have the Einstein Static Universe. This is the only solution which is 
neither expanding nor contracting. The remaining four solutions either asymptotically expand or contract 
towards or away from the Einstein Static case. Finally for A < 0, there is only one solution for which 
the universe expands and subsequently recoUapses. 

Finding the solutions when k = or —1 is similar, and the behaviour of R depends only on 
whether A is positive or negative. For A > there are two solutions, for A < 0, there is only one. When 




Fig. 1: The qualitative behaviour of the scale factor R as obtained using the condition Q < A for the three choices 
of the curvature constant k 



A = 0, we obtain the standard notion that closed universes recoUapse, while open and flat universes 
expand forever. When A 7^ 0, these simple associations are spoiled as a closed universe can expand 
forever (for large enough A) and open and flat universes can recoUapse (for A < 0). 

Exact solutions to the equations of motion can be obtained relatively easily in terms of conformal 
coordinates. We can, for example, rewrite the metric as 



ds^ = dt^ - R^{t) [dx^ + f{x) {dO^ + sin^ Odcf)'^)] 



where 



{sinhx k = —1 
X k = . 

sinx k = +1 

We can go further and define a conformal time coordinate using Rdrj = dt so that 

ds"^ = R\i]) [dt^ - dx^ + fix) {dO^ + sin^ Odf)] . 
In terms of these coordinates, the Friedmann equation becomes 

R" + kR=^{p-3p)R\ 
where ' denotes a derivative with respect to rj and is easily solved 



for p = and 



cosh r] — 1 k = —1 
R oc ( jf/2 k = 

1 — COST] k = +1 



sinhr/ k = —1 
R(x i T] k = 

sin 7/ k = +1 



t oc 



t oc 



sinh r] — r] k = —1 
rf/Q k = 

rj — sin 7/ k = +1 

cosh?7 — 1 k = —1 
7/2/2 k = 

1 — COST/ k = +1 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



for p = p/3. 

In the absence of a cosmological constant, one can define a critical energy density pc such that 
p = Pc for k = 

p^ = 3H^/8ttGn. (18) 

In terms of the present value of the Hubble parameter this is 

Pc = 1.88 X 10~2^/io^ g cm-3, (19) 

where 

ho = Ho/{WO km Mpc"^ s"^). (20) 

The cosmological density parameter is then defined by 

n = -^. (21) 

Pc 

It is useful to also define a deceleration parameter 

RqRq 



Rl 



(22) 



This standard definition was formulated under the presumption that the expansion rate of the Universe is 
in fact slowing down. As noted above, and discussed further below, modern measurements indicate the 
opposite. That is, the expansion is accelerating (meaning that go < 0). The {ii) component, Eq. (O, can 
be written in terms of go as 

-2,„H„^ = f-?^ (23) 

when the pressure is neglected. This can be combined with the Friedmann equation, Eq. ([4]), and rewritten 
as 

^ = A + F2(2go-l), (24) 

Ko 

or 

^ = Hl{^^o-q^-l). (25) 

Rq I 

Furthermore, when A = 0, go = ^o/2 so that fc = 0, +1, —1 corresponds to il = 1, J] > 1 and il < 1. 
Observational limits on Kq and fj are Ell 

/io = 0.71 ± 0.01 Oo = 1.006 ± 0.006. (26) 

It is important to note that 17 is a function of time or of the scale factor. The qualitative evolution 
of Vl is shown in Fig.[2]for A = 0. For a spatially flat Universe, Vl = \ always. When k = +1, there is a 
maximum value for the scale factor R. At early times (small values of R), Q always tends to one. Note 
that the fact that we do not yet know the sign of k, or equivalently whether i7 is larger than or smaller 
than unity, implies that we are at present still at the very left in the figure. What makes this peculiar is 
that one would normally expect the sign of k to become apparent after a Planck time of 10^^^ s. It is 
extremely puzzling that some 10^*^ Planck times later, we still do not know the sign of k. 

The Friedmann equation also lends itself to integration to determine the age of the Universe. Note 
that for a given equation of state, we can write p = po{Ro/R)^'^, and 

Hot= C ^, (27) 




Fig. 2: The evolution of the cosmological density parameter fJ as a function of the scale factor for a closed, open 
and spatially flat Universe 



where ^a = A/3Hq. When A = and Oq = 1, this is easily integrated to give 

9 

7 = 1 
4 



or 



3Ho 
1 
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(28) 



(29) 



Because of the finite age of the Universe in the Big Bang model, there is a particle horizon cor- 
responding to the maximum distance traversed by light. In general, we can write the proper distances 
between a location specified by comoving coordinate ti,ri and our location at t and r = in terms of 
the metric ([T]) 

r^ dr' 
dp = R{t) / , (30) 

^0 Vl 



For hght paths {ds^ = 0) 



* dt' 



\/l - kr'"^ 
dr' 



ii ^{'t') Jo v^l - kr 



/2 



(31) 



As ti — )■ 0, ri becomes the maximum coordinate distance from which we can receive a signal. Thus the 
particle horizon is defined by 



"^^ dr' 



dn = R{t) / -^ 
Jo Vl 



R{t) 



t j+/ 



dt' 



Vl - kr'^ ""Vo R{t') 
For A = and Oq = 1, we again obtain very simple solutions: 

dn = 'it 7 = 1 



or 



dn = 2t 7 



(32) 



(33) 



(34) 



Note that because dn/R grows with time, we see more of the Universe as time goes on. That is, new 
information is continuously entering our particle horizon. 



1.2 The hot thermal Universe 

The epoch of recombination occurs when electrons and protons form neutral hydrogen through e^ +p — )■ 
H +7 at a temperature Tr ~ few x 10'^ K ~ 1 eV. For T < Tr, photons are decoupled while for 
T > Tji, photons are in thermal equilibrium and at higher temperatures, the Universe is radiation dom- 
inated and the content of the radiation plays a very important role. Today, the content of the microwave 
background consists of photons with Tq = 2.725 it 0.001 K [3|. We can calculate the energy density of 
photons from 

p-y = E^dn^, (35) 

where the density of states is given by 

dn-y = ^[exp{E^/T) - l^^q^dq, (36) 



and g^ = 2 simply counts the number of degrees of freedom for photons, E^ = q is just the photon 
energy (momentum). (I am using units such that h = c = ks = i and will do so through the remainder 
of these lectures.) Integrating Eq. (|35] ) gives 

Py = ^T\ (37) 

which is the familiar blackbody result. In addition, we also have 



1 4p 2C(3) 3 

-P s = -— n^ = t;—' 

3^ ST ^ 7r2 



p = -p s = -— n^ = — T^T . (38) 



In general, at very early times, at very high temperatures, other particle degrees of freedom join 
the radiation background when T ~ rrij for each particle type i if that type is brought into thermal 
equilibrium through interactions. In equilibrium, the energy density of a particle type i is given by 



Eidriq^ (39) 

and 

dnq, = ^bMiEq, - Pi)/T] ± l]-\^dq, (40) 

where again gi counts the total number of degrees of freedom for type i, 

E,^ = {mj + qff\ (41) 

Pi is the chemical potential if present and ± corresponds to either Fermi or Bose statistics. 
In the limit that T ^ mi the total energy density can be conveniently expressed by 




^HE^B + ?E^H^^^-^^(^)^^ (42) 



r2 
30' 

where gB{F) are the total number of boson (fermion) degrees of freedom and the sum runs over all boson 
(fermion) states with m <^ T. The factor of 7/8 is due to the difference between the Fermi and Bose 
integrals. Equation (l42l) defines N{T) by taking into account new particle degrees of freedom as the 
temperature is raised. 

In the radiation dominated epoch, Eq. Q can be integrated (neglecting the T-dependence of N) 
giving us a relationship between the age of the Universe and its temperature 



90 



327r^GNN{T) 
6 



1/2 



-^2 



T~\ (43) 



Table 1: Effective numbers of degrees of freedom in the Standard Model 



Temperature 


New particles 


4N{T) 


T < me 


7's + i/'s 




29 


nie < T < m^ 


e± 




43 


m^ <T < rriTr 


M^ 




57 


m^ <T < T* 


vr's 




69 


Tc <T < mcharm 


- vr's + u,u. 


,d,d,s,s + gluons 


247 


rric < T < m-j- 


c,c 




289 


nir <T < mbottom 


T± 




303 


nib <T < mw,z 


6,6 




345 


mw,z <T < muiggs 


W^,Z 




381 


ruH <T < mtop 


H^ 




385 


mt<T 


t,i 




427 



*Tc corresponds to the confinement-deconfinement transition between quarks and hadrons. N{T) is shown in 
Fig.|3]for Tc = 150 and 400 MeV. 

Put into a more convenient form 

tT^,y = 2A[N{T)]~^/\ (44) 

where t is measured in seconds and TueV in units of MeV. 

The value of N{T) at any given temperature depends on the particle physics model. In the standard 
5[/(3) X SU{2) X C/(l) model, we can specify A^(r) up to temperatures of order 100 GeV. The value 
of N in the Standard Model can be seen in Table[T] 

At higher temperatures, N{T) will be model dependent. For example, in the minimal SU{5) 
model, one needs to add to N(T), 24 states for the X and Y gauge bosons, another 24 from the adjoint 
Higgs, and another 6 (in addition to the 4 already counted in W^, Z and H) from the 5 of Higgs. Hence 
for T > Mx in minimal SU{5), N{T) = 160.75. In a supersymmetric model this would at least 
double, with some changes possibly necessary in the table if the lightest supersymmetric particle has a 
mass below rrit. 

The notion of equilibrium also plays an important role in the standard Big Bang model. If, for 
example, the Universe were not expanding, then given enough time, each particle state would come 
into equilibrium with every other. Because of the expansion of the Universe, certain rates might be too 
slow indicating, for example, in a scattering process that the two incoming states might never find each 
other to bring about an interaction. Depending on their rates, certain interactions may pass in and out of 
thermal equilibrium during the course of the Universal expansion. Qualitatively, for each particle i, we 
will require that some rate Fj involving that type be larger than the expansion rate of the Universe or 

Ti > H, (45) 

in order to be in thermal equilibrium. 

A good example for a process in equilibrium at some stage and out of equilibrium at others is that 
of neutrinos. If we consider the standard neutral or charged-current interactions such as e"*" + e~ o u + u 
oxe + v-ir^e + v etc., very roughly the rates for these processes will be 

T = n{av), (46) 

where {av) is the thermally averaged weak interaction cross section 

{av)r^ O(10-2)rVM^, (47) 
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Fig. 3: The effective numbers of relativistic degrees of freedom as a function of temperature, assuming a quark- 
hadron transition temperature of 150 MeV and 400 MeV (dashed) 

and n ~ T^ is the number density of leptons. Hence the rate for these interactions is 



r„k~ o(io-2)rVM^. 



The expansion rate, on the other hand, is just 

1/2 /e^3 \l/2 



H = l'^!^') = (^iV(T)] T'^lMp ~ 1.66iV(r)V2r2/Mp. 



(48) 



(49) 



The Planck mass Mp = Gjy 



1.22 X 10^^ GeV. 
Neutrinos will be in equilibrium when F^k > H or 

T > (500M^)/Mp)i/3_ 1 MeV. 



(50) 



The temperature at which these rates are equal is commonly referred to as the decoupling or freeze-out 
temperature and is defined by 

TiTd) = H{Td). (51) 

For temperatures T > T^, neutrinos will be in equilibrium, while for T < T^ they will not. Basically, in 
terms of their interactions, the expansion rate is just too fast and they never 'see' the rest of the matter 
in the Universe (nor themselves). Their momenta will simply redshift and their effective temperature 
(the shape of their momenta distribution is not changed from that of a blackbody) will simply fall with 
T ~ 1/R. 

The relation RT ~ const is a direct consequence of the energy conservation equation Q. Indeed, 
using s = dp/dT, this equation can be rewritten as 



|(«^') = °' 



(52) 



making it more a statement of conservation of (comoving) entropy than energy (which is not conserved 
in comoving coordinates). 

Soon after neutrino decoupling, the e^ pairs in the thermal background begin to annihilate (when 
T <, rrie). Because the neutrinos are decoupled, the energy released heats up the photon background 



8 



relative to the neutiinos. The change in the photon temperature can easily be computed from entropy 
conservation. The neutrino entropy must be conserved separately from the entropy of interacting parti- 
cles. If we denote T>, the temperature of photons, and e^ before annihilation, we also have Ty = r> as 
well. The entropy density of the interacting particles at T = T> is just 



3r> 3^ 2^^30' 



5> = 0^ = 0(2 + -)(-)t3, (53) 



while at r = T<, the temperature of the photons just after e annihilation, the entropy density is 



2 



s< = -— = -(2)(— )T|, (54) 

< 3r< 3^ ^^30^ <' 

and by conservation of entropy s<^ = s> and 

(T</r>)3 = 11/4. (55) 

Thus, the photon background is at higher temperature than the neutrinos because the e^ annihilation 
energy could not be shared among the neutrinos, and 

Ty = {A/llf/^T^ ~ 1.9 K. (56) 

For further reading on standard Big Bang cosmology see Refs. |l4]-0. 

1.3 The Cosmic Microwave Background 

There has been a great deal of progress in the last several years concerning the determination of both 0.^. 
and JIa- Cosmic Microwave Background (CMB) anisotropy experiments have been able to determine 
the curvature (i.e., the sum of 0.^ and Q.j\) to better than one per cent, while observations of type la 
supernovae at high redshift and baryon acoustic oscillations provide information on (nearly) orthogonal 
combinations of the two density parameters. 

The CMB is of course deeply rooted in the development and verification of the Big Bang model 
and Big Bang Nucleosynthesis (BBN) |10|. Indeed, it was the formulation of BBN that led to the 
prediction of the microwave background lITTl . The argument is rather simple. BBN requires tempera- 
tures greater than 100 keV, which according to the Standard Model time-temperature relation, Eq. (|44]) . 
corresponds to timescales less than about 200 s. The typical cross section for the first link in the nucle- 
osynthetic chain is 

o-t;(p + 71^1) + 7) ~5 X 10"2°cmVs. (57) 

This implies that it was necessary to achieve a density 

„ IQi^ cm-^ (58) 

avt 

for nucleosynthesis to begin. The density in baryons today is known approximately from the density of 
visible matter to be hbo ~ 10^^ cm^^ and since we know that the density n scales as R^^ ^ T^, the 
temperature today must be 

To = {nBo/n)^/'^TBBN ~ 10 K (59) 

thus linking two of the most important tests of the Big Bang theory. 

Of course it was not until many years later that the microwave background radiation was discov- 
ered by Penzias and Wilson [12] while perfecting a radio antenna to track the Echo satellite. They found 
a background noise which could not be eliminated corresponding to a temperature of 3.5 it 1 K. One of 
the most important papers on modem cosmology was published with the title "A measurement of excess 
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Fig. 4: The power in the microwave background anisotropy spectrum as measured by WMAP JS), Boomerang lfT6l . 
QUaD |23|, CBI [ISJ, and ACBAR |22|. Figure courtesy of D. Scott Il24l. 



antenna temperature at 4080-Mc/s". This was followed by the seminal paper by Dicke, Peebles, Roll, and 
Wilkinson [13] putting this observation in a cosmological context. Subsequently, there have been many 
observations of the CMB culminating in the COBE observation |[T4l which determined the temperature 
to an unprecedented level, set aside any lingering doubts about the true black body nature of the CMB, 
and discovered the intrinsic anisotropics in the background. 

An enormous amount of cosmological information is encoded in the angular expansion of the 
CMB temperature 

T{e,(l)) = Y,aemYim{9,4>)- (60) 

im 

The monopole term characterizes the mean background temperature of Ty = 2.725 ± 0.001 K as deter- 
mined by COBE 1 3 1, whereas the dipole term can be associated with the Doppler shift produced by our 
peculiar motion with respect to the CMB. In contrast, the higher order multipoles are directly related to 
energy density perturbations in the early Universe. When compared with theoretical models, the higher 
order anisotropics can be used to constrain several key cosmological parameters. In the context of simple 
adiabatic cold dark matter (CDM) models, there are nine of these: the cold dark matter density, ^mh^', 
the baryon density, ^sh'^', the curvature — characterized by ritotai; the bubble parameter, h; the optical 
depth, r; the spectral indices of scalar and tensor perturbations, n^ and n^; the ratio of tensor to scalar 
perturbations, r; and the overall amplitude of fluctuations, Q. 

Microwave background anisotropy measurements have made tremendous advances in the last few 
years. The power spectrum |[T5H23l has been measured relatively accurately out to multipole moments 
corresponding to £ ~ 2000. A compilation of recent data is shown in Fig. |4] |[24]| . where the power at 



10 



a 





nil y/ y / 


/ / ' ' 


1 1 


1.0 


'i//n/ // 


/ / 


Ho 




^\ A/)/ / / 


/ 


(krTis-''Mpc-'') 


0.8 


y^/ ** // / 


^SNLS 


30 




/ / /n\ * V / / 




40 




/// '" AV / 




50 


0.6 


// ly/wi; ■ 




60 




/ /iV '/x-^ 




70 


0.4 




%.WMAP 


80 
90 




/ / m '»> 




100 


0.2 


y / 111 ' * * 


\ K"' 


i; ,•• 




/ III * * ' 


\ 




0.0 


I'/i >>>BAO \ 




1 1 1 








III * ' ' 




\ 




1 Ml Ml 


1 1 1 


\ 1 V. 



0.0 0.2 0.4 0.6 0.8 



1.0 



1.2 



1.4 



Fig. 5: Two-dimensional confidence regions in the (f^m, f^A) plane. The coloured Monte Carlo points derive from WMAP 1211 
and show that the CMB alone requires a flat Universe JIa + f^m — 1 if the Hubble constant is not too high. The SNe la 
results 1271 very nearly constrain the orthogonal combination JIa — ^m- Also shown is the region obtained from baryon 
acoustic oscillations f291. The intersection of these constraints is the most direct (but far from the only) piece of evidence 
favouring a flat model with VLm — 0.25. 



each i is given by {21 + l)Q/(47r), and Ce = (|a£m|^)- 

As indicated above, the details of this spectrum enable one to make accurate predictions of a large 
number of fundamental cosmological parameters. The results of the WMAP data (with other information 
concerning the power spectrum) is shown in Table |2] For details see Refs. 



Table 2: WMAP determinations of cosmological parameters 





WMAP alone 


WMAP + BAO + SN 


rJin/l^ 


0.133 ±0.006 


0.136 ± 0.004 


QbH^ 


0.0227 ± 0.0006 


0.0227 ± 0.0006 


h 


0.72 ±0.03 


0.705 ±0.013 


Us 


QfiS+°-°i'^ 

U.yDci_o.oi5 


0.960 ± 0.013 


T 


0.087 ±0.017 


0.084 ± 0.016 


Oa 


0.74 ±0.03 


0.726 ±0.015 



Of particular interest to us here is the CMB determination of the total density ritot as well as the 
matter density Qrn- There is strong evidence that the Universe is flat or very close to it. As noted earlier, 
the best determination of rJtotai is 1-006 ± 0.006. Furthermore, the matter density is significantly larger 
than the baryon density implying the existence of cold dark matter. Also, the baryon density, as we will 
see below, is consistent with the BBN production of D/H and its abundance in quasar absorption systems. 
The apparent discrepancy between the CMB value of r^tot and Qm, though not conclusive on its own, 
is a sign that a contribution from the vacuum energy density or cosmological constant, is also required. 
The preferred region in the Q,m — ^A plane is shown in Fig.|5] 

The presence or absence of a cosmological constant is a long standing problem in cosmology. We 
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know that the cosmological term is at most a factor of a few times larger than the cuiTent mass density. 
Thus from Eq. dJ]), we see that the dimensionless combination GnA <, 10^^^^. Also shown in Fig. [51 
are the results from SN la II26[|27I1 and baryon acoustic oscillations ESll . Taken together, we are led to a 
seemingly conclusive picture. The Universe is nearly flat with Otot — 1- However, the density in matter 
makes up only 23% of this total, with the remainder in a cosmological constant or some other form of 
dark energy. 

2 Lecture 2: Inflation and Baryogenesis 

Despite the successes of the standard Big Bang model, there are a number of unanswered questions 
that appear difficult to explain without imposing what may be called unnatural initial conditions. The 
resolution of these problems may lie in a unified theory of gauge interactions or possibly in a theory 
which includes gravity. For example, prior to the advent of grand unified theories (GUTs), the baiyon- 
to-photon ratio, could have been viewed as being embarrassingly small. Although we still do not know 
the precise mechanism for generating the baryon asymmetry, many quite acceptable models are available. 
In a similar fashion, it is hoped that a field theoretic description of inflation may resolve the problems 
outlined below. 

2.1 Cosmological problems 

2.1.1 The curvature problem 

As noted above, the determined value of Q in Eq. (I26l ) is curious since at the present time we do not know 
even the sign of the curvature term in the Friedmann equation ^, i.e., we do not know if the Universe is 
open, closed or spatially flat. 

The curvature problem (or flatness problem or age problem) can manifest itself in several ways. 
For a radiation dominated gas, the entropy density s ~ T^ and R ~ T^^. Thus assuming an adiabatically 
expanding Universe, the quantity k = k/ R^T"^ is a dimensionless constant. If we now apply the limit in 
Eq. ([261) to Eq. (El) (with A = 0) we find 



k ^ (Op - 1 )^0^ 
ii2T2 To 



k = ^^^^ = ^"" _7"° < 10-6°. (61) 



This limit on k represents an initial condition on the cosmological model. The problem then becomes 
what physical processes, if any, in the early Universe produced a value of k so extraordinarily close to 
zero (or Q. close to one). A more natural initial condition might have been k ^ 0(1). In this case the 
Universe would have become curvature dominated at T ~ lO^^Mp. For k = +1, this would signify the 
onset of recollapse. As already noted earlier, one would naturally expect the effects of curvature (seen in 
Fig. |2]by the separation of the three curves) to manifest themselves at times on the order of the Planck 
time as gravity should provide the only dimensionful scale in this era. If we view the evolution of il in 
Fig. |2]as a function of time, then it would appear that the time Iq = 13.7 Gyr = 8 x lO^'^Mp^^ (~ the 
current age of the Universe) appears at the far left of the x-axis, i.e., before the curves separate. Why 
then has the Universe lasted so long before revealing the true sign of /c? 

2.1.2 The horizon problem 

Because of the cosmological principle, all physical length scales grow as the scale factor R{t) ~ t^/^'^ , 
with 7 defined by p = (7 — l)p. However, as we have seen, there is a particle horizon dnit) ~ t as 
defined in Eq. (I32l ). For 7 > |, scales originating outside of the horizon will eventually become part of 
our observable Universe. Hence we would expect to see anisotropics on large scales f30l . 

In particular, let us consider the microwave background today. The photons we observe have 
been decoupled since recombination at T^^ ~ 3000 K. At that time, the horizon volume was simply 
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Vd oc t/, where td is the age of the Universe at T = Td- Then td = to{To/Td)^/'^ ~ 2 x 10^ yrs, where 
To = 2.725 K [3J is the present temperature of the microwave background. Our present horizon volume 
Vq oc to^ can be scaled back to td (corresponding to that part of the Universe which expanded to our 
present visible Universe) Vo(td) ^ Vo(To/T(j)'^. We can now compare Vo(td) and V^. The ratio 

Vojtd) VoTo^ tp^n^ 4 

' oc 5- OC -^ — ^ ~ 5 X 10 (62) 

Vd VdTj" tdV 

corresponds to the number of horizon volumes or casually distinct regions at decoupling which are en- 
compassed in our present visible horizon. 

In this context, it is astonishing that the microwave background appears highly isotropic on large 
scales with AT/T = l.libO.2 x 10^^ at angular separations of 10° 1 14|. The horizon problem, therefore, 
is the lack of an explanation as to why nearly 10^ causally disconnected regions at recombination all had 
the same temperature to within one part in 10^^. 

2.1.3 Density perturbations 

Although it appears that the Universe is extremely isotropic and homogeneous on very large scales (in 
fact the Standard Model assumes complete isotropy and homogeneity) it is very inhomogeneous on small 
scales. In other words, there are planets, stars, galaxies, clusters, etc. On these small scales there are large 
density perturbations namely 6p/p » 1. At the same time, we know from the isotropy of the microwave 
background that on the largest scales, 6p/p ~ 3AT/T ~ 0(10^^) |[T4l . and these perturbations must 
have grown to Sp/p ~ 1 on smaller scales. 

In an expanding Universe, density perturbations evolve with time [4]. The evolution of the Fourier 
transformed quantity -^{k,t) depends on the relative size of the wavelength A ~ /c^^ and the horizon 
scale H^^. For k <^ H, (always true at sufficiently early times) 6p/p oc t while for k ^ H, 6p/p is ~ 
constant (or grows moderately as In t) assuming a radiation dominated Universe. In a matter dominated 
Universe, on scales larger than the Jean's length scale (determined by kj = A-KGNPmatter/'Vs'^, "Vs = 
sound speed) perturbations grow with the scale factor R. Because of the growth in 6p/p, the microwave 
background limits force 6p/p to be extremely small at early times. 

Consider a perturbation with wavelength on the order of a galactic scale. Between the Planck time 
and recombination, such a perturbation would have grown by a factor of 0(10^^) and the anisotropy 
limit of 6p/p < 10^^ implies that 6p/p < 10^^^ on the scale of a galaxy at the Planck time. One should 
compare this value with that predicted from purely random (or Poisson) fluctuations of Sp/p ~ 10^^*^ 
(assuming 10®" particles (photons) in a galaxy) ||3T1 . The extent of this limit is of course related to the 
fact that the present age of the Universe is so great. 

An additional problem is related to the formation time of the perturbations. A perturbation with 
a wavelength large enough to correspond to a galaxy today must have formed with wavelength modes 
much greater than the horizon size if the perturbations are primordial, as is generally assumed. This is 
due to the fact that the wavelengths red-shift as A ~ i? ~ t^'^ while the horizon size grows linearly. It 
would appear that a mechanism for generating perturbations with acausal wavelengths is required. 

2.1.4 The magnetic monopole problem 

In addition to the much desired baryon asymmetry produced by grand unified theories, a less favourable 
aspect is also present; GUTs predict the existence of magnetic monopoles. Monopoles will be pro- 
duced |[32l whenever any simple group [such as SU{5)] is broken down to a gauge group which contains 
a U{1) factor [such as SU{3) x SU{2) x U{1)]. The mass of such a monopole would be 

Mm ~ Mgut/oigut ~ 10^^ GeV. (63) 



13 



The basic reason monopoles are produced is that in the breaking of SU(5), the Higgs adjoint needed 
to break SU(5) cannot align itself over all space [33]. On scales larger than the horizon, for example, 
there is no reason to expect the direction of the Higgs field to be aligned. Because of this randomness, 
topological knots are expected to occur and these are the magnetic monopoles. We can then estimate that 
the minimum number of monopoles produced [34] would be roughly one per horizon volume or causally 
connected region at the time of the SU{5) phase transition tc, 

rim ~ (2tc)"^ (64) 

resulting in a monopole-to-photon ratio expressed in terms of the transition temperature of 

rijn (lQTc\ .... 

(65) 



r^ 



The overall mass density of the Universe can be used to place a constraint on the density of monopoles. 
For Mjn ~ 10^*^ GeV and O^/iq^ ;^ 1 we have that 

^< 0(10-25). (66) 

rij 



The predicted density, however, from Eq. (1651 ) for Tc ~ Mqut is 

^ ~ 10-9. (67) 

Hence we see that standard GUTs and cosmology have a monopole problem. 

2.2 Inflation 

All of the problems discussed above can be neatly resolved if the Universe underwent a period of cos- 
mological inflation 0511361 . That is, if the Universe at some stage becomes dominated by the vacuum as 
could be the case during a phase transition, our assumptions of an adiabatically expanding Universe may 
not be valid. Indeed, we expect several cosmological phase transitions to have occurred including the 
breakdown of a Grand Unified symmetry such as SU(5) — ;■ SU(3)cX SU(2)j;,x U(l)y or the electroweak 
transition SU(2)lx U(l)y — )• U(l)em. or possibly some other non-gauged transition. 

During a phase transition, the motion of a scalar field will be described by a scalar potential. If 
the solution to the equations of motion for the scalar field leads to a slowly evolving scalar field (this 
will depend on the details of the potential), the Universe may become dominated by the vacuum energy 
density associated with the potential near the initial field value, say w 0. The energy density of the 
symmetric vacuum V{^) acts as a cosmological constant with 

A = 87ry(0)/Mp2. (68) 

During this period of slow evolution, the energy density due to radiation or matter will fall below the 
vacuum energy density, p <^ V{0). When this happens, the expansion rate will be dominated by the 
constant V{0) and from Eq. ^ we find the exponentially expanding solution given in Eq. ^. When 
the field evolves towards the global minimum it will begin to oscillate about the minimum, energy will 
be released during its decay and a hot thermal Universe will be restored. If released fast enough, it will 
produce radiation at a temperature Tr'^ <, V{0). In this reheating process, entropy has been created and 
{RT)f > {RT)i. Thus we see that during a phase transition, the relation RT ^ constant need not hold 
true and our dimensionless constant k may actually not have been constant. 

If during the phase transition the value of RT changed by a factor of C(10^*^), the cosmological 
problems discussed above would be solved. The isotropy would in a sense be generated by the immense 
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expansion; one small causal region would get blown up and our entire visible Universe would have 
been at one time in thermal contact. In addition, the parameter k could have started out 0(1) and 
have been driven small by the expansion. The wavelengths of density perturbations would have been 
stretched by the expansion A ~ i? making it appear that A S> H^^ or that the perturbations have left the 
horizon. Rather, it is the size of the causally connected region that is no longer simply H^^. However, 
not only does inflation offer an explanation for large scale perturbations, it also offers a source for the 
perturbations themselves ll37l . Monopoles would also be diluted away. 

The cosmological problems could be solved if 

Ht > 65 (69) 

where r is the duration of the phase transition. In a successful theory, density perturbations are produced 
and do not exceed the limits imposed by the microwave background anisotropy, the vacuum energy 
density was converted to radiation so that the reheated temperature is sufficiently high, and baryogenesis 
is realized. 

In the original (old) inflationary scenario, the phase transition determined by a potential with a 
large barrier separating the false and true vacua proceeds via the formation of bubbles [38]. The Universe 
reheats with the release of entropy which must occur through bubble collisions and the transition is 
completed when the bubbles fill up all of space. It is known |[39]| . however, that the requirement for 
a long timescale r is not compatible with the completion of the phase transition. The Universe as a 
whole remains trapped in the exponentially expanding phase containing only a few isolated bubbles of 
the broken SU{3) x SU{2) x f/(l) phase. 

The well-known solution to the dilemma of old inflation is called the new inflationary scenario BOl . 
New inflation (it was at the time) was originally based on symmetry breaking using a flat potential of the 
Coleman-Weinberg form Ref. ||4TI . Instead of proceeding by tunnelling and the formation of bubbles, 
the transition takes place more or less uniformly on large scales. The details of the inflationary transition 
are determined from the equations of motion. 

A Lagrangian for a scalar field which includes a scalar potential may be incorporated into the total 
action including gravity 



/ = / d^x^ { £2 - ^d^'i^d'^cl) - VicP) ] . (70) 



E_ 1 

2^~ 2 

The equation of motion for a scalar field (p can be derived from the energy-momentum tensor 



T^u = df,^d,<p - \g^,,dp<i)dP<i) - g^,,V{(t>). (71) 
By associating p = Tqo and p = R~'^{t)Tii we have 

/> = ^02 ^ii?-2(t)(V0)2 + !/(</.), (72) 

p = \4>^- \R-\t){V(pf - F(0), (73) 

2, D 

and from Eq. ^ we can write the equation of motion (by considering a homogeneous region, we can 
ignore the gradient terms) 

ip + ?,H4) = -dVjdcp. (74) 

Consider now the approximation dV jdcp ~ {d'^V / d4>'^)4>; the equation of motion becomes H2l 

<?i + 3iJ0 + m2(0)(/) = O, (75) 
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with m?{cl)) = d'^V/dcjP' < 0. The solution when \m?\ » H^ grows exponentially as ~ e'™'* while 
for \m?\ <^ H^ the scalar field grows as ^ ~ e'™' */3^ ju (jjg latter case the field moves very slowly 
during a time period 

r~3i//|mp. (76) 

This approximation is known as the slow-rollover approximation. 

If the scalar mass is tuned somewhat, m^ ~ 10^ GeV, a significant amount of inflation is possible. 
From Eq. (1761 ) one sees that 



H' V 



4 



i4 



Ht ~ ~ ~ 10^ (11) 

for v ~ 10^^ GeV. Reheating no longer occurs via the collisions of bubbles, but by the decay of scalar 

field oscillations. As the scalar field settles to its minimum, the solution to the equations of motion look 

like 

(pit) r^ —sin mt (78) 

mt 

and the reheat temperature is 

Tr ~ {TDM^f^ Vd < Hj, (79) 

where F/) is the scalar field decay rate and Hj is the value of H during inflation. 

In addition to producing il = 1, which is clearly seen from Eq. (|6TI ) as /c — ;■ 0, new inflation is 
capable of producing scale invariant density perturbations [37| of the type preferred for galaxy formation 
models. However, the original BOl new inflationary models based on a Coleman-Weinberg 11411 type 
of SU{5) breaking produced density fluctuations with magnitude dp/p ~ O(IO^) rather than 6p/p ~ 
10^^ as needed to remain consistent with microwave background anisotropics. Other more technical 
problems |43| concerning slow rollover and the effects of quantum fluctuations also pass doom on this 
original model. 

General models of inflation can be described by a few so-called slow roll parameters e and r/. 
These are given by 



^ + 1^ = ±l(±WMI(YM. 



and 



For sufficient inflation both of these parameters must be small. The amount of inflation (given by the 
total number of e-foldings, A^) is given by e, 

As noted above, inflation leads to a nearly scale-free spectrum of density fluctuations, A(fc) oc 
6p/p with a power spectrum of the form A^(A;) oc fc"^^. The spectral index is determined by the 
inflationary potential 

n ~ 1 - 6e + 2r/. (83) 

In addition to the scalar perturbations, tensor perturbations (gravity waves) will also be produced during 
inflation with spectral index Ug ~ — 2e. The ratio of the amplitudes of the tensor to scalar perturbations 
is an important observable also given in terms of the slow roll parameter, e, r ~ We. 

A generic model of inflation can be described by a potential of the form 

V{ct)) = p^P{ct)), (84) 
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Table 3: Slow roll parameters for chaotic inflationary models 



e 


1/120 


1/60 


V 


1/120 


1/40 


n 


0.97 


0.95 


r 


0.13 


0.27 



where </> is the scalar field driving inflation, the inflaton, n is an as yet unspecified mass parameter, 
and P(0) is a function of which possesses the features necessary for inflation, but contains no small 
parameters. That is, P{(l)) takes the form 

P{ri) = P(0) + rn^rf + X^rf + X^rj'^ + ..., (85) 

where all of the couplings in P are 0(1) and ... refers to possible non-renormalizable terms. Most of the 
useful inflationary potentials can be put into the form of Eq. (l84l ). 

The requirements for successful inflation boil down to: 1) enough inflation; and 2) density pertur- 
bations of the right magnitude. The latter reduces approximately to 

^~O(100)-^. (86) 

For large scale fluctuations of the type measured by COBE 1. 14.1 . we can use Eq. (I86l ) to fix the inflationary 
scale /I of the inflaton potential B4l : 



/^^ 



^ = few X 10"*. (87) 

Fixing {p? /Mp) has immediate general consequences for inflation f45l. For example, the Hubble 
parameter during inflation, H^ ~ (87r/3)(/x'^/Mp) so that H r^ lO^^Mp. The duration of inflation is 
r ~ Mp/fi'^, and the number of e-foldings of expansion is Ht ~ SiT{Mp/ ji'^) ~ 10^. If the inflaton 
decay rate goes as F ~ m^/Mp ~ fi^/Mp, the Universe recovers at a temperature Tr ~ (FAfp)^/^ ~ 
p^/M^ ~ 10-i^Mp ~ 10^ GeV. 

Two commonly studied potentials are associated with a form of inflation known as chaotic infla- 
tion B6i . In these models it is assumed that as part of an initially chaotic state > Mp with V{cl)) ~ Mp. 
Once these assumptions are made, chaotic models of inflation are by far the simplest. Typical models for 
chaotic inflation in terms of a single scalar field are described by the following scalar potentials 1146 [1471 



V = ^Xcl)^ (88) 

or 

V = ]-m^(t)^. (89) 

That's all! Nothing more complicated is necessary. It is assumed that at the Planck time, all fields (f){x) 
satisfy V{ct)) < Afp^ and {d^cf))^ < Mp'^. It is also assumed that there exist domains sufficiently large 
I > H^^ with (f) homogeneous and cj) > Mp. 

For suitably large initial values of the scalar field 0, the Universe expands quasiexponentially. 
Sufficient inflation requires only that (po > few xMp. Although this is not a strong constraint on chaotic 
models, a stronger constraint is derivable from the consideration of density perturbations. One finds that 
6p/p < 10~^ for A < 10~^^ or m < lO^^Mp. The slow roll parameters are easily determined for these 
two models and are given in Table [3l 

Finally, CMB measurements can be used to test inflationary models by determining or limiting the 
slow roll parameters. For example, WMAP ||2| is able to set a constraint in the r, n plane as shown. For 
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example, atn = 0.95, the 68% (95%) CL upper limit on r is < 0.07 (0.17), while at n = 0.97, r < 0.18 
(0.28). As one can see, while the m'^cp'^ model is well within the constraints, the A^^ model is not. 

2.3 Baryogenesis 

It appears that there is apparently very little antimatter in the Universe. To date, the only antimatter 
observed is the result of a high-energy collision, either in an accelerator or in a cosmic-ray collision 
in the atmosphere. There has been no sign to date of any primary antimatter, such as an anti-helium 
nucleus a found in cosmic -rays. In addition, the number of photons greatly exceeds the number of 
baryons. Indeed, the value of il^^^ as determined by WMAP 121 and listed in Table |2] corresponds to a 
baryon-to-photon ratio of 

ry = !^£^I^ ^ !^ :^ 27AnBh^ = (6.23 ± 0.17) x IQ-^o. (90) 

In the Standard Model, the entropy density today is related to rij by 

27r^ 21 4 

s = rT(2 + )n^ = 7.04n^, (91) 

90C(3)^ 4 11^ ^ ^' 

so that Eq. (l90l ) implies ub/s r^ 8.8 x 10^^^. In the absence of baryon number violation or entropy pro- 
duction this ratio is conserved, however, and hence represents a potentially undesirable initial condition. 

Let us for the moment assume that in fact r] = 0. We can compute the final number density of nu- 
cleons left over after annihilations of baryons and antibaryons have frozen out. At very high temperatures 
(neglecting a quark-hadron transition) T > 1 GeV, nucleons were in thermal equilibrium with the photon 
background and ub = Ub = (3/2)n^ (a factor of 2 accounts for neutrons and protons and the factor 3/4 
for the difference between Fermi and Bose statistics). As the temperature fell below m]\j annihilations 
kept the nucleon density at its equilibrium value {uB/n^) = (vr^/^(mAr/T)^/^/2'^/^C(3))exp(— rriAr/T) 
until the annihilation rate Ta — nB'm^'^ fell below the expansion rate. This occurred at T ~ 20 MeV. 
However, at this time the nucleon number density had already dropped to 

ub/ti^ = ng/n^ ~ 10"^^, (92) 

which is eight orders of magnitude too small Il48l aside from the problem of having to separate the 
baryons from the antibaryons. If any separation did occur at higher temperatures (so that annihilations 
were as yet incomplete) the maximum distance scale on which separation could occur is the causal 
scale related to the age of the Universe at that time. At T =20 MeV, the age of the Universe was 
only t = 2 X 10^^ s. At that time, a causal region (with distance scale defined by 2ct) could only 
have contained W~^Mq which is very far from the galactic mass scales, lO^^Af0, we are asking for 
separations to occur. 

2.3.1 The out-of-equilibrium decay scenario 

The production of a net baryon asymmetry requires baryon number violating interactions, C and CP 
violation and a departure from thermal equilibrium B9l . The first two of these ingredients are contained 
in GUTs, the third can be realized in an expanding Universe where, as we have seen, it is not uncommon 
that interactions come in and out of equilibrium. In SU(5), the fact that quarks and leptons are in the 
same multiplets allows for baryon non-conserving interactions such as e~ +(i -H- n + u, etc., or decays of 
the supermassive gauge bosons X and Y such as X — )• e~ + d, ii + u. Although today these interactions 
are very ineffective because of the very large masses of the X and Y bosons, in the early Universe when 
T ~ Mx ~ 10^^ GeV these types of interactions should have been very important. C and CP violation is 
very model dependent. In the minimal SU(5) model, as we will see, the magnitude of C and CP violation 
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is too small to yield a useful value of rj. The C and CP violation in general comes from the interference 
between tree level and first loop corrections. 

The departure from equilibrium is very common in the early Universe when interaction rates 
cannot keep up with the expansion rate. In fact, the simplest (and most useful) scenario for baryon 
production makes use of the fact that a single decay rate goes out of equilibrium. It is commonly referred 
to as the out-of-equilibrium decay scenario |50|. The basic idea is that the gauge bosons X and Y (or 
Higgs bosons) may have a lifetime long enough to insure that the inverse decays have already ceased so 
that the baryon number is produced by their free decays. 

More specifically, let us call X, either the gauge boson or Higgs boson which produces the baryon 
asymmetry through decays. Let a be its coupling to fermions. For X a gauge boson, a will be the GUT 
fine structure constant, while for X a Higgs boson, (47ra)^'^ will be the Yukawa coupling to fermions. 
The decay rate for X will be 

Yd ^ aMx. (93) 

However, decays can only begin occurring when the age of the Universe is longer than the X lifetime 
r^^, i.e., when Td > H 

aMx > N{Tfl^T'^/Mp, (94) 

or at a temperature 

T^ < aMxMpN{T)~^/^. (95) 

Scatterings, on the other hand, proceed at a rate Ts ~ a^T^/M"^ and hence are not effective at lower 
temperatures. To be in equilibrium, decays must have been effective as T fell below Mx in order to 
track the equilibrium density of X's (and X's). Therefore, the out-of-equilibrium condition is that at 
T = Mx,Td <Hor 

Mx ^ aMp{N{Mx)y^'^ ~ lO^^a GeV. (96) 

In this case, we would expect a maximal net baryon asymmetry to be produced. 

To see the role of C and CP violation, consider the two channels for the decay of an X gauge 
boson: X — )• (l)uu, {2)e^ d. Suppose that the branching ratio into the first channel with baryon number 
B = —2/3 is r and that of the second channel with baryon number B = +1/3 is 1 — r. Suppose in 
addition that the branching ratio for X into {\)uu with baryon number B = +2/3 is f and into {2)e^d 
with baryon number B = —1/3 is 1 — f. Though the total decay rates of X and X (normalized to 
unity) are equal as required by CPT invariance, the differences in the individual branching ratios signify 
a violation of C and CP conservation. 

Denote the parity (P) of the states (1) and (2) by t or |, then we have the following transformation 
properties: 

Under CPT: r(X ^ 1 1) = T{l_i^ X) 
Under CP : T{X ^ 1 t) = T{X -^ 1 i) (97) 

Under C: r(X ^ 1 t) = T{X^l\). 

We can now denote 

r = T{X ^ 1 t) + r(X ^ 1 i) (98) 

f = r(x^it) + r(x^ 11). (99) 

The total baryon number produced by an X, X decay is then 

2 1, ,2 1, 

A5 = r + -(l-r) + -f (1-r) 

= f-r = r(x^ It) + r(x^ 11) -r(x^ It) -r(x^ 11). (loo) 

One sees clearly therefore, that from Eqs. ( |97l) if either C or CP are good symmetries, IS.B = 0. 
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Fig. 6: The time evolution of the baryon asymmetry with B = L = initially 



In the out-of-equilibrium decay scenario 11501 . the total baryon asymmetry produced is proportional 
to AB = {f — r). If decays occur out of equilibrium, then at the time of decay nx ~ n^ at T < Mx- 
We then have 

nB_ ^ {AB)nx _ {AB)nx _ ,^_2, 
s s 



10-'^{AB). 



(101) 



s N{T)n^ 

The schematic view presented above can be extended to a complete calculation given a specific model ISTl 
l52l . see also Ref. Il53l for reviews. 

The time evolution for the generation of a baryon asymmetry is shown in Fig.[6l As one can see, for 
large values of Mx, i.e., values which satisfy the lower limit given in Eq. (I96l ). the maximal value for the 
baryon asymmetry ub/ s ~ 10~^e is achieved. This confirms numerically the original out-of-equilibrium 
decay scenario ISOll . For smaller values of Mx an asymmetry is still produced which, however, is smaller 
due to partial equilibrium maintained by inverse decays. The growth of the asymmetry as a function of 
time is now damped, and it reaches its final value when inverse decays freeze out. Finally, by studying 
different initial conditions, one can show that the result for the final baryon asymmetry is in fact largely 
independent of the initial baryon asymmetry. 

From Eq. (IIOII ) it is clear that a complete calculation of n^/s will require a calculation of the CP 
violation in the decays (summed over parities) which we can parametrize by 



T{X ^ 1) - T{X -^ 1) Imr 
T{X ^l) + T{X ^ 1) ^ Ref 



(102) 



At the tree level, as one can see, T{X — ;• 1) oc gl^g^ is real and there is no C or CP violation. At the one 
loop level for gauge boson decay, there is also no net contribution to e, and we must turn to Higgs decay. 
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Fig. 7: One-loop contribution to the C and CP violation with two Higgs five-plets 



At least two Higgs five-plets are required to generate sufficient C and CP violation. With two five-plets, 
H and H', the interference of diagrams of the type in Fig. |7] will yield a non- vanishing e 154)1 . 

e oc Im(a'"^a6'6l') / , (103) 

if the couplings a ^ o! and 6/6'. 

There are of course many alternative methods to generate the baryon asymmetry, though each 
makes use of the same three ingredients. For example, a supersymmetric mechanism proposed by Affleck 
and Dine |[55l makes use of flat directions in the scalar potential. There are many such flat directions, 
and some of these yield an non-vanishing expectation value to GUT baryon number violating operators. 
Supersymmetry breaking perturbs the flatness, leading to the cosmological evolution of the scalar which 
oscillates about the global (charge and colour conserving) minimum of the potential. Baryon number is 
stored in these oscillations and a net asymmetry is produced as these fields decay. 

Another mechanism to generate the baryon asymmetry employs the heavy right-handed neutrinos 
used in the see-saw mechanism to generate neutrino masses f56l. The simplest of such mechanisms is 
based on the decay of a right-handed neutrino-like state [57]. This mechanism is certainly novel in that 
it does not require grand unification at all. By simply adding to the Lagrangian a Dirac and Majorana 
mass term for a new right-handed neutrino state, 

£ 9 Mv^v" + XHLv", (104) 

the out-of-equilibrium decays z^"^ — )■ L + H* and v^ ^ L* + H will generate a non-zero lepton number 
L / 0. The out-out-equilibrium condition for these decays translates to 10~^A^Mp < M and M could 
be as low as 0(10) TeV. (Note that once again in order to have a non-vanishing contribution to the C 
and CP violation in this process at 1-loop, at least two flavours of z/^ are required. For the generation of 
masses of all three neutrino flavors, three flavours of u'^ are required.) Electroweak sphaleron effects ||58]| 
can transfer this lepton asymmetry into a baryon asymmetry. If sphalerons are in equihbrium, the baryon 
number can be expressed in terms of B — L 

OS 

B = -{B-L). (105) 

In the absence of a primordial B — L asymmetry, the baryon number is erased by equilibrium processes. 
Right-handed neutrinos produce a net lepton asymmetry and hence a net B — L yielding the final baryon 
asymmetry given by Eq. (11051 ). 

3 Lecture 3: Big Bang Nucleosynthesis 

The standard model |[59l of Big Bang Nucleosynthesis (BBN) is based on the relatively simple idea of 
including an extended nuclear network into a homogeneous and isotropic cosmology. Apart from the 
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input nuclear cross sections, the theory contains only a single parameter, namely the baryon-to-photon 
ratio, 7], and even that has been fixed by WMAP |2|. The theory then allows one to make predictions 
(with well-defined uncertainties) of the abundances of the light elements, D, ^He, ''He, and ^Li |[60ll . 

Conditions for the synthesis of the light elements were attained in the early Universe at tempera- 
tures T >, I MeV. In the early Universe, the energy density was dominated by radiation with 

P=^(^ + l + l^-^'^^^ (106) 

from the contributions of photons, electrons and positrons, and N^ neutrino flavours (at higher tem- 
peratures, other particle degrees of freedom should be included as well). At these temperatures, weak 
interaction rates were in equilibrium. In particular, the processes 

n + e"*" -f-)- p + I'e 
n + Vf, <-)• p + e~ 

n o p + e~ +Ue (107) 

fix the ratio of number densities of neutrons to protons. At T » 1 MeV, [n/p] ~ 1. 

As we have seen in the case of neutrino interactions, the weak interactions do not remain in 
equilibrium at lower temperatures. Freeze-out occurs when the weak interaction rate T^k ~ G\T^ falls 
below the expansion rate which is given by the Hubble parameter H ~ T^/Mp. The /3-interactions 
in Eq. (11071) freeze out at about 0.8 MeV. As the temperature falls and approaches the point where the 
weak interaction rates are no longer fast enough to maintain equilibrium, the neutron-to-proton ratio is 
given approximately by the Boltzmann factor, {n/p) ~ g-^rn/T ^ -^1^^ where Am is the neutron- 
proton mass difference. After freeze-out, free neutron decays drop the ratio slightly to about 1/7 before 
nucleosynthesis begins. A useful semi-analytic description of freeze-out has been given 116 111621 . 

The nucleosynthesis chain begins with the formation of deuterium by the process, p + n — )■ D + 7. 
However, because of the large number of photons relative to nucleons, 77^^ = n^/uB ~ 10^'', deuterium 
production is delayed past the point where the temperature has fallen below the deuterium binding energy, 
Eb = 2.2 MeV (the average photon energy in a blackbody is Er^ ~ 2.7 T). This is because there are 
many photons in the exponential tail of the photon energy distribution with energies E > Eb despite the 
fact that the temperature or E^ is less than Eb- The degree to which deuterium production is delayed can 
be found by comparing the qualitative expressions for the deuterium production and destruction rates, 

Tp « UBav (108) 

When the quantity r]~^ex.p{—EB/T) ~ 1, the rate for deuterium destruction (D + 7 — )• p + n) finally 
falls below the deuterium production rate and the nuclear chain begins at a temperature T ~ 0.1 MeV. 

The dominant product of Big Bang nucleosynthesis is ^He and its abundance is very sensitive to 
the (n/p) ratio 

2 (n/p) 
[1 + (n/p)] 

i.e., an abundance of close to 25% by mass. Lesser amounts of the other light elements are produced: 
D and ^He at the level of about 10^^ by number, and ^Li at the level of 10"^*^ by number. The gap at 
A = 8 prevents the production of other isotopes in any significant quantity. The nuclear chain is shown 
in Fig. [i 

Historically, BBN as a theory explaining the observed element abundances was nearly abandoned 
due its inability to explain all element abundances. Subsequently, stellar nucleosynthesis became the 
leading theory for element production ||63L However, two key questions persisted. 1) The abundance of 

22 



Yp = ^^ ; /7^^ w 0.25, (109) 




Fig. 8: The nuclear network used in BBN calculations 



^He as a function of metallicity is nearly flat and no abundances are observed to be below about 23% as 
exaggerated in Fig. |9l In particular, even in systems in which an element such as oxygen which traces 
stellar activity is observed at extremely low values (compared with the solar value of 0/H « 4.9 x 10^^), 
the ^He abundance is nearly constant. This is very different from all other element abundances (with 
the exception of ^Li as we shall see below). For example, in Fig. \T0\ the N/H vs. 0/H correlation is 
shown |[64l . As one can clearly see, the abundance of N/H goes to 0, as 0/H goes to 0, indicating a 
stellar source for nitrogen. 2) Stellar sources cannot produce the observed abundance of D/H. Indeed, 
stars destroy deuterium and no astrophysical site is known for the production of significant amounts of 
deuterium Il65l . Thus we are led back to BBN for the origins of D, ^He, ^He, and ^Li. 



3.1 Abundance predictions 

Because standard BBN theory rests upon the Standard Model of particle physics, the electroweak aspects 
of the calculation are very well-determined and do not introduce an appreciable uncertainty. Instead, the 
major uncertainties come from the thermonuclear reaction rates. There are 11 key strong rates (as well 
as the neutron hfetime) which dominate the uncertainty budget Il66l - l69l . In contrast to the situation for 
much of stellar nucleosynthesis, BBN occurs at high enough temperatures that laboratory data exist at 
and even below the relevant energies, so that no extrapolation is needed. Monte Carlo techniques II66II67I 
are used to determine the best-fit abundances, and their uncertainties, at each rj. 

Recently the input nuclear data have been carefully reassessed 1166116814711 . leading to improved 
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Fig. 9: The ''He mass fraction as determined in extragalactic H II regions as a function of O/H 
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Fig. 10: The nitrogen and oxygen abundances in the same extragalactic H II regions with observed ""^He shown in 
Fig.S 



precision in the abundance predictions. In addition, polynomial fits to the predicted abundances and 
the error correlation matrix have been given II72II731 . The NACRE Collaboration presented an updated 
nuclear compilation ITTOI . For example, notable improvements include a reduction in the uncertainty in 
the rate for ^Re{n,p)T from 10% ^ to 3.5% and for T(a, jYU from ~ 23-30% CI to ~ 4%. Since 
then, new data and techniques have become available, motivating new compilations. Within the last year, 
several new BBN compilations have been presented I 73ir75]477l . 

The light element abundances are shown in Fig. [TT]as a function of r] ITTTII . The plot shows the 
abundance of ^He by mass Y and the abundances of the other three isotopes by number. The curves 
indicate the central predictions from BBN, while the bands correspond to the uncertainty in the predicted 
abundances. The uncertainty range in ^He reflects primarily the la uncertainty in the neutron lifetime. 

In standard BBN with N^ = 3, the only free parameter is the density of baryons (strictly speaking, 
nucleons), which sets the rates of the strong reactions. Because standard BBN is a one-parameter theory, 
any abundance measurement determines rj, while additional measurements overconstrain the theory and 
thereby provide a consistency check. BBN has thus historically been the premier means of determining 
the cosmic baryon density. 
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Fig. 11: The predictions of standard BBN [771 showing the primordial abundances as a function of the baryon-to- 
photon ratio 77. Abundances are quantified as ratios to hydrogen, except for ^Hc which is given in baryonic mass 
fraction Yp = pne/pB- The lines give the mean values, and the surrounding bands give the la uncertainties. 



The release of the first-year WMAP results on the anisotropy spectrum of the CMB were a land- 
mark event for all of cosmology, but particularly for BBN. As discussed above, the value of r/ has been 
fixed by CMB measurements as given by Eq. (|90l ). Thus, within the context of the Standard Model, BBN 
becomes a zero-parameter theory, and the light element predictions are completely determined to within 
the uncertainties in rj and the BBN theoretical errors. Comparison with light element observations then 
can be used to restate the test of BBN-CMB consistency, or to turn the problem around and test the 
astrophysics of post-BBN light element evolution iTTSl . 



3.2 Light element observations and comparison with theory 

BBN theory predicts the abundances of D, •^He, ^He, and ^Li, which are essentially determined at 
t ~ 180 s. Abundances are, however, observed at much later epochs, after stellar nucleosynthesis 
has commenced. The ejecta from stellar processing can alter the light element abundances from their 
primordial values, but also produce heavy elements such as C, N, O, and Fe ('metals'). Thus one seeks 
as trophy sical sites with low metal abundances, in order to measure light element abundances which are 
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closer to primordial. For all of the light elements, systematic errors are an important and often dominant 
limitation to the precision of the primordial abundances. 

3.2.1 D/H 

In recent years, high-resolution spectra have revealed the presence of D in high-redshift, low-metallicity 
quasar absorption systems (QAS), via its isotope-shifted Lyman-a absorption. These are the first mea- 
surements of light element abundances at cosmological distances. It is believed that there are no astro- 
physical sources of deuterium ll65l . so any measurement of D/H provides a lower limit to primordial D/H 
and thus an upper limit on rj. Recent observations by FUSE show a wide dispersion in the deuterium 
abundance in local gas seen via its absorption, (D/H)iocai gas = (0.5 — 2.2) x 10~^ |791. This surpris- 
ingly large spread, taken together with the positive correlation of D/H with temperature and metallicity 
along various sightlines, led |[79l to suggestions that deuterium may suffer significant and preferential 
depletion onto dust grains. In this case the true local interstellar D/H value would he at the upper limit 
of the observed values, giving (D/H)jgj^ ^ (2.31 it 0.24) x 10~^. However, extracting a primordial 
deuterium value requires a Galactic chemical evolution model (e.g., Ref. [80]), whose model depen- 
dences yield uncertainties in the determination of the primordial deuterium abundance. Many of these 
models do not predict significant D/H depletion JSTll at high redshift, and in this case the high-redshift 
measurements are expected to recover the primordial deuterium abundance. 

The deuterium abundance at low metallicity has been measured in several quasar absorption sys- 
tems 1,82,1 . The weighted mean value of the seven systems with reliable abundance determinations is 
log D/H = —4.55 ± 0.03 where the error includes a scale factor of 1.72 and corresponds to D/H = 
(2.82 lb 0.21) X 10^^. These are shown in Fig.[l2l Since the D/H shows considerable scatter it is likely 
that systematic errors dominate the uncertainties. In this case it may be more appropriate to derive the 
uncertainty using sample variance (see, for example, Ref. |66|) which gives a more conservative range 
log D/H = -4.55 ± 0.08 or D/H = (2.82 ± 0.53) x 10"^. 

Using the WMAP value for the baryon density (l90l) . the primordial D/H abundance is predicted to 
be 1771 

(D/H)p = (2.49 ± 0.17) X 10~^ (110) 

As one can see from Fig. \TT\ this is in good agreement with the average of the seven best determined 
quasar absorption system abundances noted above, particularly when systematic uncertainties are taken 
into account. 

3.2.2 ^He 

We observe ^He in clouds of ionized hydrogen (H II regions), the most metal-poor of which are in dwarf 
galaxies. There is now a large body of data on ^He and CNO in these systems [ 83 ] . These data confirm 
that the small stellar contribution to helium is positively correlated with metal production. Recently a 
careful study of the systematic uncertainties in '^He, particularly the role of underlying absorption, has 
led to a higher value for the primordial abundance of ^He |[84l . Using a subset of the highest quality 
from the data of Izotov and Thuan f83l, all of the physical parameters listed above including the ^He 
abundance were determined self-consistently with Monte Carlo methods 1.85.1 . The extrapolated ^He 
abundance was determined to be Y^ = 0.249 it 0.009. Conservatively, it would be difficult at this time 
to exclude any value of Yp inside the range 0.232-0.258. 

At the WMAP value for r], the '^He abundance is predicted to be 11771 

Yp = 0.2486 ± 0.0002. (Ill) 

This is in excellent agreement with the most recent analysis of the "^He abundance [84]. Note also that the 
large uncertainty ascribed to this value indicates that while ^He is certainly consistent with the WMAP 
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Fig. 12: D/H abundances shown as a function of [Si/H]. Labels denote the background quasi-stellar objects (QSO), 
except for the local interstellar value (LISM; 1791 ). 



determination of the baryon density, it does not provide for a highly discriminatory test of the theory at 
this time. 



3.2.3 ^Li 

The systems best suited for Li observations are metal-poor stars in our Galaxy. Observations have long 
shown f86l that Li does not vary significantly in Pop II stars with metallicities < 1/30 of solar — the 
'Spite plateau'. Precision data suggest a small but significant correlation between Li and Fe [87] which 
can be understood as the result of Li production from Galactic cosmic rays [88]. Extrapolating to zero 
metallicity one arrives at a primordial value [[89l Li/H|p = (1.23lo;5g) x 10"^°. 

Figure [13] shows the different Li components for a model with (^Li/H)p = 1.23 x 10"^*^ as a 
function of the iron abundance expressed as the log of Fe/H relative to the solar value. The linear slope 
produced by the model is independent of the input primordial value. The model [90J includes, in addition 
to primordial ^Li, lithium produced in Galactic cosmic-ray nucleosynthesis (primarily a + a fusion), 
and ^Li produced by the z^-process during type II supemovae. As one can see, these processes are not 
sufficient to reproduce the population I abundance of ^Li (at near solar [Fe/H] <, 0), and additional 
production sources are needed. 
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Fig. 13: Contributions to the total predicted lithium abundance from the adopted Galactic chemical evolution 
model of Ref. |90|, compared with low metallicity stars and a sample of high metallicity stars. The solid curve is 
the sum of all components. 



A recent reanalysis of the '^He(a, 7)^Be reaction, which is the most important ^Li production 
process in BBN, was considered in detail in Ref. [91]. When the new rate is used a high ^Li abundance 
is found liZZl at the WMAP value of r/ 
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(112) 



^Li/H = (5.24^ 

This represents a 23% increase in ^Li over previous calculations iT/Bl . The increase is primarily due to an 
increase in the ■^He(a, 7)^Be cross section. Newer data fOT] implies 17% increase in this reaction leading 
to a 16% increase in ^Li. In addition, the 1.5% increase in j] from the 3-year to 5-year WMAP data [21 
leads to a 3% increase in ^Li and finally another 1% increase is due to updated pn rates. In addition, the 
uncertainty in the BBN ^Li abundance is roughly a factor of 2 times smaller than previous determinations. 
This value for primordial ^Li is in clear contradiction with most estimates of the primordial Li abundance. 
Several attempts at explaining this discrepancy by adjusting some of the key nuclear rates have proved 
unsuccessful Il7ni92ll9^. 

An important source for systematic error lies in the derived effective temperature of the star. [Li] 
= log(^Li/H) + 12 is very sensitive to the temperature, with 9[Li]/5Teff ~ 0.065 - 0.08. Unfortunately 
there is no standard for determining effective temperatures, and, for a given star, there is considerable 
range depending on the method used. This spread in temperatures was made manifest in the recent work 
of Melendez and Ramirez 1941 using the infra-red flux method (IRFM) which showed differences for 
very low metallicities ([Fe/H] < -3) by as much as 500 K, with typical differences of ~ 200 K with 
respect to that of Ref. [871. As a consequence the derived ^Li abundance was significantly higher with 
Li/H|p = (2.34 ± 0.32) x 10-i° f94,95]. 

Recently a dedicated set of observations were performed with the specific goal of determining the 
effective temperature in metal-poor stars ll96l . Using a large set of Fe I excitation lines (~ 100 lines 
per star), the Boltzmann equation was used with the excitation energies, Xi to determine the temperature 
through the distribution of excited levels. Again, there was no evidence for the high temperatures reported 
in Ref. [94], rather, temperatures were found to be consistent with previous determinations. The mean 
^Li abundance found in Ref. [|96l was Li/H = (1.3 - 1.4 ± 0.2) x 10"^°, consistent with the bulk of 
prior abundance determinations. 
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There are of course other possible sources of systematic uncertainty in the ^Li abundance. It 
is possible that some of the surface ^Li has been depleted if the outer layers of the stars have been 
transported deep enough into the interior, and/or mixed with material from the hot interior; this may 
occur due to convection, rotational mixing, or diffusion. Estimates for possible depletion factors are 
in the range r^ 0.2-0.4 dex, i.e., factors of 1.6-2.5 [97|. Recent attempts to deplete the ^Li abundance 
through diffusion introduce a source of turbulence tuned to fit the abundances of heavy elements in 
NGC6397 ||98]| . Once parameters are set, the degree of lithium depletion becomes a prediction of the 
model. For this cluster, a depletion factor of 0.25 dex is found, i.e., a factor of 1.8. Note that while this 
depletion factor would bring the previous BBN result of ^Li/H = 4.26 x 10"^'^ [i73il to a value close the 
value observed in that cluster [99.1 (2.2 x 10^^'^), a larger depletion factor is needed with the new BBN 
value for ^Li given above. It is also not clear whether this mechanism will work for the wide range of 
stellar parameters seen in the field. As noted above, the Li data show a negligible intrinsic spread in 
Li. Any mechanism which reduces significantly the abundance of ^Li must do so uniformly over a wide 
range of stellar parameters (temperature, surface gravity, metallicity, rotational velocity, etc.). 

It is also possible that the lithium discrepancy is a sign of new physics beyond the Standard Model. 
One possibility is the cosmological variation of the fine structure constant. Varying a would induce a 
variation in the deuterium binding energy and could yield a decrease in the predicted abundance of 
^Li [100|. A potential solution to both lithium problems is particle decay after BBN which could lower 
the ^Li abundance (and produce some ^Li as well) | I101| . This has been investigated in the framework 
of the constrained minimal supersymmetric Standard Model if the lightest supersymmetric particle is 
assumed to be the gravitino II 1021 and indeed, some models have been found which accomplish these 



goals IT031 . 

3.2.4 ^He 

Since ^He is also quite sensitive to the baryon density, one might hope that it too could be used as 
a baryometer. Observations of H II regions in our own Galaxy yield values of the '^He/H ratio that 
are compatible with calculations of the primordial value II1041I1051 . However, the extrapolation from 
the observations to a primordial abundance is complicated by the unknown chemical evolution of ^He. 
Indeed, one does not even know whether '^He/H is increasing or decreasing with cosmic time. Thus a 



primordial extrapolation yields only an order-of-magnitude range of allowable values of '^He/H II106II . 

3.3 Beyond the Standard Model 

Given the simple physics underlying BBN, it is remarkable that it still provides one of the most effective 
tests for physics beyond the Standard Model. Limits on particle physics beyond the Standard Model come 
mainly from the observational bounds on the ^He abundance. As discussed earlier, the neutron-to-proton 
ratio is fixed by its equilibrium value at the freeze-out of the weak interaction rates at a temperature 
Tj ~ 1 MeV modulo the occasional free neutron decay. Furthermore, freeze-out is determined by the 
competition between the weak interaction rates and the expansion rate of the Universe 



Gf%' ~ r„eak(T/) = HiTf) ~ ^/GNNTf^. (113) 

In the Standard Model, the number of relativistic particle species at 1 MeV is A^ = ^ + jN,y. The 
presence of additional neutrino flavours (or any other relativistic species) at the time of nucleosynthesis 
increases the overall energy density of the Universe and hence the expansion rate leading to a larger 
value of Tf, (n/p), and ultimately Yp. Because of the form of Eq. (11131 ) it is clear that just as one can 
place limits 111071 on N,y, any changes in the weak or gravitational coupling constants can be similarly 
constrained (for a discussion see Ref. BIOS II ). 

The helium curve in Fig. \TT\ was computed taking Ni, = 3; the computed abundance scales as 
^Ybbn — 0.013AA^^ Ii6ril . The dependence of the light element abundances on N^, is shown in 
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Fig. 14: BBN abundance predictions fTSl as a function of the baryon-to-photon ratio rj, for N,y = 2-7. The bands 
show the la error bars. Note that for the isotopes other than Li, the error bands are comparable in width to the 
thickness of the abundance curve shown. All bands are centred on N^ = 3. 



Fig. [H EH- For a fixed value of r/ = (6.14 ± 0.25) x 10"^° (slightly below the current WMAP 
value) and Yp = 0.249 ± 0.009, the likelihood distribution for N^ is shown by the shaded region in 
Fig. [T5]||109i Also shown for comparison are the likelihood distribution based the WMAP value of t] us- 
ing D/H alone, Yp and D/H, and the result based on BBN alone. Despite the increased uncertainty in the 
He abundance, it still provides the strongest constraint on N^. D/H is nonetheless becoming competitive 
in its ability to set limits on N^. 

The 95 % CL upper limits to N,^ are given in Tabled In all cases the preferred values for N,^ are 
consistent with N^, and in many cases are much closer to N,^ than la. This restates the overall consis- 
tency among standard BBN theory, D and '^He observations, and CMB anisotropics. It also constrains 
departures from this scenario. The combined limit using BBN + light elements -i- CMB limit is II 10911 



2.67 < iVv < 3.85 



(114) 



at 68% CL. 

4 Lecture 4: Dark Matter 

Evidence for dark matter in the Universe is available from a wide range of observational data. As dis- 
cussed many times above, the analysis of the cosmic microwave background anisotropics leads to the 
conclusion that the curvature of the Universe is close to zero indicating that the sum of the fractions 
of critical density, Q, in matter and a cosmological constant (or dark energy) is very close to one [2]. 
When combined with a variety of data including results from the analysis of type la supernovae observa- 
tions I26II27II and baryon acoustic oscillations |[28l one is led to the concordance model where $!„ ~ 0.23 
and r^A ~ 0.73 with the remainder (leading to O^tot = 1) in baryonic matter. This is in addition to the 
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Fig. 15: The likelihood disttibution for N^ based on the WMAP value of ?; and Yp (shaded), WMAP and D/H 
(dashed), WMAP and both Yp and D/H^i (dotted). Also shown is the result without imposing the WMAP value for 
7y (long dashed). 

Table 4: This table shows constraints placed on N^ and rj by various combinations of observations. Shown are 
the 68% confidence limits determined by marginalizing the 2-D likelihood distribution. Also shown are the 95% 
upper limits on SN^, = N^, — 3, given that SN^, > 0. 



Observations rjio = lO^^rj 


N^ 


J^ i>,max 
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classic evidence from galactic rotation curves IIIIOII . which indicate that nearly all spiral galaxies are 
embedded in a large galactic halo of dark matter leading to rather constant rotational velocities at large 
distances from the centre of the galaxy (in contrast to the expected f ^ ~ 1/r behaviour in the absence of 
dark matter). Other dramatic pieces of evidence can be found in combinations of X-ray observations and 
weak lensing showing the superposition of dark matter (from lensing) and ordinary matter from X-ray 
gas milH and from the separation of baryonic and dark matter after the collision of two galaxies as seen 
in the Bullet cluster II112II . For a more complete discussion see Ref. II113L 

From the first column of Table|2j we can obtain the density of cold dark matter from the difference 
between the total matter density and the baryon density [2] 
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0.1099 ±0.0062 



(115) 



or a 2cj range of 0.0975-0.1223 for VtcDuh'^- 
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Fig. 16: Summary plot 111 171 of the relic density of Dirac neutrinos (solid) including a possible neutrino asymmetry 
of ?7,, =: 5 X 10"" (dotted) 



4.1 Neutrinos 

Dark matter must be both long-lived or stable and electrically and colour neutral. As a result, once 
baryons and neutrinos are eliminated as candidates, one must look beyond the Standard Model. From 
WMAP, we already know that the baryon density is far below the requisite amount in cold dark matter. 
Light neutrinos (m < 30 eV) are a long-time standard when it comes to non-baryonic dark matter HI 141 . 
Light neutrinos are, however, ruled out as a dominant form of dark matter because they produce too much 
large scale structure II115II . The energy of density of light neutrinos with my < 1 MeV can be expressed 
at late times as p^ = jirriyn^. Imposing the constraint Vlyh^ <, 0.12, translates into a strong constraint 
(upper bound) on Majorana neutrino masses 1111611 : 



m-tot 



^m,, < 11 eV, 



(116) 



where the sum runs over neutrino mass eigenstates. The Umit for Dirac neutrinos depends on the in- 
teractions of the right-handed states. The limit (11161 ) and the corresponding initial rise in Q.yh? as a 
function of my is displayed in Fig. [161 Much stronger limits on the sum of neutrino masses are possible 
when combining the WMAP data with large scale structure surveys. A typical limit is TTitot < 0.7 eV or 
Vtyh^ < 0.0076 |2|. 

The calculation of the relic density for neutrinos more massive than ~ 1 MeV, is substantially 
more involved. The relic density is now determined by the freeze-out of neutrino annihilations which 
occur at T < m,^, after annihilations have begun to seriously reduce their number density [118|. For 
particles which annihilate through approximate weak scale interactions, annihilations freeze out when 

T ~ m^/20. 

Based on the leptonic and invisible width of the Z boson, experiments at LEP have determined that 
the number of neutrinos is N^ = 2.994 it 0.012 [119]. Thus, LEP excludes additional neutrinos (with 
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standard weak interactions) with masses rrii^ < 45 GeV. The mass density of ordinary heavy neutrinos 
is found to be very small, f^j^/i^ < 0.001 for masses m„ > 45 GeV up to m,y ~ 0(100) TeV 11181. 
Laboratory constraints for Dirac neutrinos are available II120L excluding neutrinos with masses between 
10 GeV and 4.7 TeV. This is significant, since it precludes the possibility of neutrino dark matter based 
on an asymmetry between u and u Iil21il . 

4.2 Axions 

Owing to space limitations, the discussion of axions as a dark matter candidate will be very brief. Ax- 
ions are pseudo-Goldstone bosons which arise in solving the strong CP problem Ill22[|123ll via a global 
U(l) Peccei-Quinn symmetry. The invisible axion II123I is associated with the flat direction of the spon- 
taneously broken PQ symmetry. Because the PQ symmetry is also explicitly broken (the CP violating 
OFF coupling is not PQ invariant), the axion picks up a small mass similar to a pion picking up a mass 
when chiral symmetry is broken. We can expect that rua ~ rn-j^f.,^/ fa where fa, the axion decay con- 
stant, is the vacuum expectation value of the PQ current and can be taken to be quite large. If we write 
the axion field as a = faO, near the minimum, the potential produced by QCD instanton effects looks 
like V ~ m'^O'^fa- The axion equations of motion lead to a relatively stable oscillating solution. The 
energy density stored in the oscillations exceeds the critical density II 12411 unless fa ^ 10^^ GeV. 

Axions may also be emitted from stars and supernovae II125I . In supernovae, axions are produced 
via nucleon-nucleon bremsstrahlung with a coupling g^N oc rriN / fa- As was noted above, the cosmo- 
logical density limit requires fa ^ 10^^ GeV. Axion emission from red giants imply 1.126,1 fa ^ 10^*^ GeV 



(though this limit depends on an adjustable axion-electron coupling), the supernova limit requires II 12711 
/a >, 2 X 10^^ GeV for naive quark model couplings of the axion to nucleons. Thus only a narrow 
window exists for the axion as a viable dark matter candidate. 

4.3 Neutralinos 

Supersymmetry is one of the best-motivated proposals for physics beyond the Standard Model. It is 
well known that supersymmetry could help stabilize the mass scale of electroweak symmetry break- 
ing by cancelling the quadratic divergences in the radiative corrections to the mass-squared of the 
Higgs boson [1281. In addition, including supersymmetric partners of Standard Model particles in the 
renormalization-group equations (RGBs) for the gauge couplings of the Standard Model would permit 
them to unify 111291 . whereas unification would not occur if only the Standard Model particles were 
included in the RGBs. 

To construct the supersymmetric Standard Model f lBOl we start with the complete set of chiral 
fermions needed in the Standard Model, and add a scalar superpartner to each Weyl fermion so that each 
field in the Standard Model corresponds to a chiral multiplet. Similarly we must add a gaugino for each of 
the gauge bosons in the Standard Model making up the gauge multiplets. The Minimal Supersymmetric 
Standard Model (MSSM) [131] is defined by its minimal field content (which accounts for the known 
Standard Model fields) and minimal superpotential necessary to account for the known Yukawa mass 
terms. As such we define the MSSM by the superpotential 

W = tij {yeii{Ve' + vaRIQU' + y..i/^Q%^] + eijixR\Rl . (117) 

In Bq. (I117I ). the indices {ij} are SU(2)2, doublet indices. The Yukawa couplings y are all 3 x 3 matrices 
in generation space. Note that there is no generation index for the Higgs multiplets. Colour and genera- 
tion indices have been suppressed in the above expression. There are two Higgs doublets in the MSSM. 
This is a necessary addition to the Standard Model which can be seen as arising from the holomorphic 
property of the superpotential. That is, there would be no way to account for all of the Yukawa terms for 
both up-type and down-type multiplets with a single Higgs doublet. To avoid a massless Higgs state, a 
mixing term must be added to the superpotential. 
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In defining tlie MSSM, we have limited the model to contain a minimal field content: the only 
new fields are those which are required by supersymmetry. Consequently, apart from superpartners, only 
the Higgs sector was enlarged from one doublet to two. Moreover, in writing the superpotential (11171 ). 
we have also made a minimal choice regarding interactions. We have limited the types of interactions 
to include only the minimal set required in the Standard Model and its supersymmetric generalization. 
There are, however, additional superpotential terms which are consistent with gauge invaiiance. These 
would lead to rapid baryon and/or lepton number violation and can be eliminated by imposing a discrete 
symmetry on the theory called i2-parity 111321 . This can be represented as 



R 



-1 



:iB+L+2s 



(118) 



where B, L, and s are the baryon number, lepton number, and spin, respectively. It is easy to see that, 
with the definition (11181 ). all the known Standard Model particles have i?-parity +1. For example, the 
electron has B = 0, L = —1, and s = 1/2, and the photon has B = L = and s = 1, so in both cases 
R = I. Similarly, it is clear that all superpartners of the known Standard Model particles have R = —1, 
since they must have the same value of B and L as their conventional partners, but differ by 1/2 unit of 
spin. If i?-parity is exactly conserved, then the additional superpotential terms must be absent from the 
theory. An immediate result of imposing i?-parity is the stability of the lightest R = —I sparticle making 
it a potential dark matter candidate. Possible choices for the lightest supersymmetric particle (LSP) are 
the neutralino, sneutrino, and gravitino. Here, I will focus only on the former. 

There are four neutralinos, each of which is a linear combination of the R = —I neutral fermions I1133II : 
the wino W^, the partner of the third component of the SU{2)i gauge boson; the bino, B; and the two 
neutral Higgsinos, Hi and if 2- The mass and composition of the LSP are determined by the gaugino 
masses, fi, and tan /?. In general, neutralinos can be expressed as a linear combination 



X = aB + PW'^ + -fHi + 6H2. 



(119) 



The solution for the coefficients a, j3, 7 and 5 for neutralinos that make up the LSP can be found by 
diagonalizing the mass matrix 
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where Mi(M2) is a soft supersymmetry breaking term giving mass to the U(l) (SU(2)) gaugino(s). 

The relic density of neutralinos depends on additional parameters in the MSSM beyond Mi , M2 , fi, 
and tan/3. These include the sfermion masses rrir and the Higgs pseudo-scalar mass tua- To determine 
the relic density it is necessary to obtain the general annihilation cross-section for neutralinos. In much of 
the parameter space of interest, the LSP is a bino and the annihilation proceeds mainly through sfermion 
exchange. 

In its generality, the minimal supersymmetric standard model (MSSM) has over 100 undetermined 
parameters. There are good arguments based on grand unification [:129.l and supergravity [134] which 
lead to a strong reduction in the number of parameters. I will assume several unification conditions 
placed on the supersymmetric parameters. In all models considered, the gaugino masses are assumed to 
be unified at the GUT scale with value mi/2 ^s are the trilinear couplings with value Aq. Also common 
to all models considered here is the unification of all soft scalar masses set equal to ttiq at the GUT scale. 
With this set of boundary conditions at the GUT scale, we can use the radiative electroweak symmetry 
breaking conditions by specifying the ratio of the two Higgs vacuum expectation values, tan /3, and the 
mass Mz to predict the values of the Higgs mixing mass parameter /i and Higgs pseudoscalar mass 
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Fig. 17: RG evolution of the mass parameters in the CMSSM. I thank Toby Falk for providing this figure. 



niA- The sign of fi remains free. This class of models is often referred to as the constrained MSSM 
(CMSSM) 11135141391 In the CMSSM, the solutions for fj. generally lead to a hghtest neutralino which 
is very nearly a pure B. 

In Fig.[T7J an example of the renormalization group running of the mass parameters in the CMSSM 
is shown. Here, we have chosen m^ /2 = 250 GeV, mo = 100 GeV, tan P = 3, Aq = 0, and /i < 0. 
Indeed, it is rather amazing that, from so few input parameters, all of the masses of the supersymmetric 
particles can be determined. The characteristic features that one sees in the figure are, for example, 
that the coloured sparticles are typically the heaviest in the spectrum. This is due to the large positive 
correction to the masses due to 03 in the RGBs. Also, one finds that the B is typically the lightest 
sparticle. But most importantly, notice that one of the Higgs mass^ goes negative triggering electroweak 
symmetry breaking II 1401 . (The negative sign in the figure refers to the sign of the mass squared, even 
though it is the mass of the sparticles which is depicted.) 

For a given value of tan (3, Aq, and sgn{fi), the resulting regions of parameter space with accept- 
able relic density and which satisfy the phenomenological constraints can be displayed on the 1111/2 — mo 
plane. In Fig. [TSla). the light shaded region corresponds to that portion of the CMSSM plane with 
tan /3 = 10, ^0 = 0, and /i > such that the computed relic density yields the WMAP value given 
in Eq. (11151 ) fT38l. The bulk region at relatively low values of mi/2 and ttiq, tapers off as m,i/2 is in- 
creased. At higher values of ttiq, annihilation cross sections are too small to maintain an acceptable relic 
density and ri^/i^ is too large. Although sfermion masses are also enhanced at large mi/2 (due to RGB 
running), co-annihilation processes between the LSP and the next lightest sparticle (in this case the f) 
enhance the annihilation cross section and reduce the relic density. This occurs when the LSP and NLSP 
are nearly degenerate in mass. The dark shaded region has m^ < m^ and is excluded. The effect of 
co-annihilations is to create an allowed band about 25-50 GeV wide in mo for mi/2 ^ 950 GeV, or 
1711/2 ~ 400 GeV, which tracks above the nif-^ = m^ contour II141L 

Also shown in Fig. [TSl a) are the relevant phenomenological constraints. These include the LEP 
limits on the chargino mass: m,^± > 104 GeV Iil42 l: on the selectron mass: mg > 99 GeV II143I : 
and on the Higgs mass: nih > 114 GeV 111441 . The former two constrain mi/2 and nio directly via 
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Fig. 18: The {mi/2,mQ) planes for (a) ta.n(3 = 10 and /i > 0, assuming ^o = 0,mt = 175 GeV and 
mh(?7ifc)|^j = 4.25 GeV. The near-vertical (red) dot-dashed lines are the contours rrih — 114 GeV, and the 
near-vertical (black) dashed line is the contour m^± — 104 GeV. Also shown by the dot-dashed curve in the lower 
left is the corner excluded by the LEP bound of wg > 99 GeV. The medium (dark green) shaded region is excluded 
by 6 — >^ S7, and the light (turquoise) shaded area is the cosmologically preferred region. In the dark (brick red) 
shaded region, the LSP is the charged fi. The region allowed by the E821 measurement of a^t at the l-cr level, is 
shaded (pink) and bounded by solid black lines, with dashed lines indicating the l-a ranges. In (b), tan /3 = 50. 



the sparticle masses, and the latter indirectly via the sensitivity of radiative corrections to the Higgs 
mass to the sparticle masses, principally m^^^. FeynHiggs |145| is used for the calculation of nih- 
The Higgs limit imposes important constraints principally on m^ /2 particularly at low tan /3. Another 
constraint is the requirement that the branching ratio for 6 — )• S7 be consistent with the experimental 
measurements [.146.1 . These measurements agree with the Standard Model, and therefore provide bounds 
on MSSM particles II147II such as the chargino and charged Higgs masses, in particular. Typically, the 
6 — )• 57 constraint is more important for fi < 0, but it is also relevant for ^u > 0, particularly when 
tan/3 is large. The constraint imposed by measurements of 6 — )• 57 also excludes small values of mi/2. 
Finally, there are regions of the (7711/2, "^o) plane that are favoured by the BNL measurement [.148ll of 
g^ — 2 at the 2-a level, corresponding to a deviation from the Standard Model calculation II149II . 

Another mechanism for extending the allowed regions in the CMSSM to large rriy^ is rapid anni- 
hilation via a direct-channel pole when m^ ~ ^tua II1351I137I . Since the heavy scalar and pseudoscalar 
Higgs masses decrease as tan /3 increases, eventually 2m^ ~ niA yielding a 'funnel' extending to large 
1711/2 and uiQ at large tan /?, as seen in Fig. fTSl b). As one can see, the impact of the Higgs mass constraint 
is reduced (relative to the case with tan /3 = 10) while that of 6 — ;■ S7 is enhanced. 

Shown in Fig. [19] are the WMAP Unes II138II of the {'mi/2,mo) plane for ^u > and values of 
tan/3 from 5 to 55, in steps A(tan/3) = 5. We notice immediately that the strips are considerably 
narrower than the spacing between them, though any intermediate point in the (?ni/2> ^n-o) plane would 
be compatible with some intermediate value of tan /3. The right (left) ends of the strips correspond to 
the maximal (minimal) allowed values of ?tii/2 and hence m-^. The lower bounds on 1111/2 are due to the 
Higgs mass constraint for tan /3 < 23, but are determined by the 5 — )• 57 constraint for higher values of 
tan/3. 

Finally, there is one additional region of acceptable relic density known as the focus-point re- 
gion [,150,1 . which is found at very high values of mo- An example showing this region is found in 
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Fig. 19: The strips display the regions of the {mi/2, mo) plane that are compatible with the WMAP determination 
of il-^h^ and the laboratory constraints for fi > and tan /3 = 5,10,15, 20, 25, 30, 35, 40, 45, 50, 55. The parts of 
the strips compatible with g^ — 2 at the 2-a level have darker shading. 



Fig- ESI plotted for tan/3 = 10, ^ > 0, and mt = 172.4 GeV. As mo is increased, the solution for fi at 
low energies as determined by the electroweak symmetry breaking conditions eventually begins to drop. 
When fi < 1711/2, the composition of the LSP gains a strong Higgsino component and as such the relic 
density begins to drop precipitously. As ttiq is increased further, there are no longer any solutions for fi. 
This occurs in the shaded region in the upper left corner of Fig.[20l The position of the focus point strip 
is very sensitive to the value of mt II151I . 

As seen in Fig. [TH the relic density constraint is compatible with relatively large values of mi/2 
and m,Q. However, all values of m.1/2 and mo are not equally viable when the available phenomenological 
and cosmological constraints are taken into account. A global likelihood analysis enables one to pin down 
the available parameter space in the CMSSM. One can avoid the dependence on priors by performing a 
pure likelihood analysis as in Ref. II152II . or a purely x^-based fit as done in Refs. II153II1541 . Here we 
present results from one such analysis II155II1561 . which used a Markov-Chain Monte Carlo (MCMC) 
technique to explore efficiently the likelihood function in the parameter space of the CMSSM. A full list 
of the observables and the values assumed for them in this global analysis are given in Ref. [, 154,1 . and 
updated in Refs. II155II156I . 

The 68% and 95% confidence-level (CL) regions in the (r)^i/2 7 w-o) plane of the CMSSM is shown 
in Fig. |2T][155]. Also shown for comparison are the physics reaches of ATLAS and CMS with 1/fb of 
integrated luminosity Iil57lll58l . (MET stands for missing transverse energy, SS stands for same-sign 
dilepton pairs, and the sensitivity for finding the lightest Higgs boson in cascade decays of supersym- 
metric particles is calculated for 2/fb of data.) The likelihood analysis assumed /i > 0, as motivated by 
the sign of the apparent discrepancy in g^ — 2, but sampled all values of tan /? and Ao'. the experimental 
sensitivities were estimated assuming tan /3 = 10 and Ao = 0, but are probably not very sensitive to 
these assumptions. The global maximum of the likelihood function (indicated by the black dot) is at 
mi/2 = 310 GeV, mo = 60 GeV, Ao = 240 GeV, tan/3 = 11, and x /Ndof = 20.4/19 (37% probabil- 
ity). It is encouraging that the best-fit point lies well within the LHC discovery range, as do the 68% and 
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Fig. 20: As in Fig. [18] showing the rT^i/2, n^o plane extended to high values of mo 



most of the 95% CL regions. 

Improvements in sensitivity have made it possible for direct detection experiments Ill59[|160ll to 
be competitive with other phenomenological constraints. The elastic cross section for x scattering on a 
nucleus can be decomposed into a scalar (spin-independent) and a spin-dependent part. Each of these 
can be written in terms of the cross sections for elastic scattering off individual nucleons. The scalar part 
of the cross section can be written as 



fsi 



Ami 



vr 



[ZU + {A-Z)f^f 



where m^ is the ^-nuclear reduced mass and 
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(122) 



for A^ = p or n. The parameters /j, are defined by 



ruNf^''^ = {N\mqqq\N) = m,Sf ), 



(123) 



and the 03^ contain the individual quark-neutralino scattering cross sections, see Refs. 111611 - 116311 for 
further details regarding the calculation of the cross section. 

The elastic scattering of neutralinos on nucleons is very sensitive to the stiangeness contribution 
to the nucleon mass and can be characterized by the parameter y which is also related to the vr-nucleon 
sigma term S^jv by 

y = S — r-5- = 1 ~ cro/S,riv , (124) 

where o"o is the change in the nucleon mass due to non-zero u and d masses and is estimated from octet 
baryon mass differences to be cjo = 36 MeV II164II . and the latest determination of T,.,^,^ = 64 MeV. 
The effects of varying these assumptions are discussed in the context of the CMSSM in Refs. 111621 
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11631 . Lattice calculations are now reaching the stage where they may also provide useful information on 
EttAt 111651 . and a recent analysis would suggest a lower value T^^^n < 40 II1661 . 

In Fig. l22l' a) we show CMSSM spin-independent neutralino-nucleon cross section, as obtained 
in a scan over all CMSSM parameters with 5 < tan/3 < 55, 100 < mi/2 < 2000 GeV, GeV < 
mo < 2000 GeV, and — 3m,i/2 < ^o !^ 3mi/2 11671 . We also allow both positive and negative /x, 
except for large tan /3 > 30, where convergence becomes difficult in the ^u < case. At low m^ < 
300 GeV, cross sections generally exceed 10^^ pb, and the largest scalar cross sections, which occur for 
m^ ~ 100 GeV, are already excluded by CDMS II |T59l and/or XENONIO f\Mi . These exclusions 
occur primarily in the focus-point region at large tan /3. On the other hand, for m^ > 400 GeV, scalar 
cross sections are well below 10~^ pb, and come from the co-annihilation strip or the rapid-annihilation 
funnel that appears at large tan /3 in the CMSSM. The effective cross sections shown are suppressed for 
points with J7^ <^ ^cdm, and there may be cancellations at larger m^ that suppress the cross sections 
substantially. These regions of parameter space will not be probed by direct detection experiments in the 
near future [ 168, 169|. The corresponding 68% and 95% CL regions in the cross section-neutralino mass 
plane from the frequentist analysis of Ref. 115611 are shown in Fig.l22lb). 
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